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Abstract 

We analyze the low-lying spectrum and eigenmodes of lattice Dirac 
operators with a twisted mass term. The twist term expels the eigen¬ 
values from a strip in the complex plane and all eigenmodes obtain a 
non-vanishing matrix element with 75 . For a twisted Ginsparg-Wilson 
operator the spectrum is located on two arcs in the complex plane. 
Modes due to non-trivial topological charge of the underlying gauge 
field have their eigenvalues at the edges of these arcs and obey a rem¬ 
nant index theorem. For configurations in the confined phase we find 
that the twist mainly affects the zero modes, while the bulk of the 
spectrum is essentially unchanged. 
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Lattice fermions with a twisted mass term have gained a lot of attention 
recently niiiEiEi. Particularly appealing for numerical simulations is the 
fact that a twisted mass term cures the notorious problem with exceptional 
configurations at a very low cost. This allows for simulations with relatively 
small pion mass close to the physical point |^. 

For small pion masses the role of topological configurations of the gauge 
field and the corresponding dynamics of low-lying modes of the Dirac opera¬ 
tor become more important for QCD phenomenology. Of particular interest 
is the index theorem |S] and its manifestation on the lattice. While for chi¬ 
ral lattice Dirac operators, obeying the Ginsparg-Wilson equation [^, an 
exact index theorem holds [Zj, for other formulations only an approximate 
realization of the index theorem can be expected. For Wilson and staggered 
fermions this has been studied numerically BU- For twisted mass fermions 
the interplay of topology and low-lying Dirac eigenmodes has not yet been 
analyzed. In this letter we close this gap, combining analytical arguments 
and numerical results. 

The twisted mass Dirac operator V for two mass-degenerate flavors of fer¬ 
mions, u and d, has the form 

V = {Do + m)l 2 + 75 T 3 , To = diag(l,-l), 

with I 2 and T 3 acting in flavor space. Alternatively one can define the Dirac 
operators for u and d separately by = D{fi) and = D{—fi) with 

D{n) = Do + m + in75 . (1) 

Do is a standard lattice Dirac operator and for the numerical results pre¬ 
sented here is chosen to be Wilson’s Dirac operator, n is the twisted mass 
parameter and m the regular mass. In our notation we assume ^ > 0. In 
this work the mass parameter m is adjusted such that for /r = 0 massless 
quarks are described (“maximal twist”). We remark, however, that for the 
study of eigenvalues and eigenvectors m is an irrelevant parameter, since it 
only shifts the whole spectrum in the complex plane. 

As remarked, QCD with a twisted mass term is particularly interesting 
for lattice simulations since exceptional configurations are avoided because 

det[D„] det[Dd] = det[D(/i)] det[D(-/r)] = det[ 75 D(/i) 75 ] det[D(-/i)] = 
det[(Dj-|-m-|-i^ 75 )(Do-|-m-i// 75 )] = dQi[{Do+m)\Do+m) + , ( 2 ) 

showing that the determinant for two flavors is strictly positive. For stag¬ 
gered fermions the regular mass term already excludes exceptional config¬ 
urations. For Wilson fermions fluctuations of real eigenvalues can become 
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large and only for relatively high values of the regular mass are exceptional 
configurations ruled out. When using the Wilson operator as input in the 
overlap projection, these fluctuations can also make the evaluation of the 
square root problematic. 

In the second line of Eq. m we have used the 75-hermiticity of the Dirac 
operator Dq, 

75 Dq 75 = Dq . ( 3 ) 

In passing we note that for non-vanishing twist parameter // this equation 
is replaced by the generalized 75-hermiticity relation m 

( 4 ) 

Since the similarity transformation on the l.h.s. leaves the eigenvalues in¬ 
variant, we conclude from (0 that the spectrum of D(—/x) is the complex 
conjugate of the spectrum of D(/x). This property is reflected in the fact, 
that the product of determinants in 0 is real. 

Let us begin our analysis of the twisted mass Dirac operator D(^) by com¬ 
paring the spectrum with and without twist. In Fig.^we show the low-lying 
eigenvalues for /x = 0 (l.h.s. plot) and fi = 0.02 (r.h.s.) on the same quenched 
configuration, generated with the Liischer-Weisz action m on a 16 ^ lattice 
at /3 = 8 . 45 . The lattice spacing is a = 0.094 fm as determined from the 
Sommer parameter in m- The eigenvalues A were computed with the im¬ 
plicitly restarted Arnold! method m- In the plots the numbers next to 
the symbols representing the eigenvalues A are the chiralities of the corre¬ 
sponding eigenvectors tpx, i.e., the matrix elements {'ipx, ysV'A)- The plots are 
representative for a larger ensemble of configurations we studied at several 
values of /x, ranging from /x = 0.01 to /U = 0.1. 

We first discuss the spectrum for zero twist (l.h.s. plot). The configu¬ 
ration we consider has topological charge Q = +1 and thus we find a real 
eigenvalue near the origin. This “topological eigenvalue” is represented by 
a filled circle, while for the other eigenvalues we use open diamonds. For a 
lattice Dirac operator obeying the Ginsparg-Wilson equation the topologi¬ 
cal eigenvalue would be exactly at the origin. For Wilson’s Dirac operator 
75-hermiticity implies that real eigenvalues near the origin take over the role 
of the zero-modes, since only for real eigenvalues r the corresponding eigen¬ 
vector has non-trivial chirality, (V’rjTsV’r) 7^ 0 . This is evident from the 75 
matrix elements displayed next to the symbols for the eigenvalues. 

Upon turning on the twist parameter the spectrum changes considerably 
(r.h.s. plot). The eigenvalues are expelled from a strip along the real axis 
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Figure 1 : Spectrum of the Wilson Dirac operator in the complex plane for 
vanishing twist (l.h.s. plot) and /i = 0.02 (r.h.s.). For both plots the same 
quenched charge Q = +1 configuration was used. The filled circle is used 
for the topological mode (zero-mode), while open diamonds represent bulk 
modes. The numbers next to the lowest eigenvalues are the values of the 75 
matrix elements (chirality) for the corresponding eigenvectors. 


excluding eigenvalues A with |ImA| < /r. Also the topological mode is shifted 
outside the “forbidden strip.” It is shifted downwards since, according to the 
index theorem, it has negative chirality for our Q = -|-1 configuration. The 
eigenvalues of topological modes with positive chirality are shifted upwards. 
This demonstrates that for /r ^ 0 the spectrum is not symmetric with respect 
to reflection at the real axis. We stress that this holds not only for the 
topological eigenvalue, but also the other eigenvalues do not come in complex 
conjugate pairs. For // = 0 the spectrum is symmetric as a consequence of 
®. Another drastic change is the fact that for fi ^ 0 all eigenvectors acquire 
non-vanishing chirality. They have positive chirality for eigenvalues in the 
upper half and negative chirality in the lower half (for positive fi). The 
absolute value of the 75 matrix element decreases with increasing distance 
from the real axis. 

Let us now try to understand these observations analytically. The 75 matrix 
element of an eigenvector V’a can be transformed as 

A(V’a,75V’a) = (V’a,75-D(m)V'a) = (V'a,-D(-^)^75V’a) = (-D(-^)V’a,75V’a) = 
(P(M)-*2/i75]V’A,75V'A) = A*(V’A,75V’A)+i2//(V’A,V'A) = A*(V'a, 75V'a)+*2//. 
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In the second step we have used (| 1 J and in the last step that the eigenvectors 
are normalized to 1 . From this equation follows 

(^/’a, 75 ' 0 a) = ( 5 ) 

establishing that all eigenvectors have non-vanishing chirality which de¬ 
creases monotonically with ImA. Using the fact that 75 is bounded, i.e., 
|(V’a,75V’a)| < 1 , we find (/r > 0) 

I ImA I > /r . (6) 

This result for the imaginary part is of course equivalent to the lower bound 
for the determinant in Eq. ©• 

We remark that the fluctuations of the real part of the eigenvalues are 
not reduced by the twisted mass term. This is obvious from comparing the 
scatter of the eigenvalues in horizontal direction for the two plots in Fig. ^ 
and was observed also on the larger sample of configurations we analyzed. 
For /i = 0 the range of these fluctuations is related to the Aoki phase. The 
fact that the horizontal spread of the eigenvalues is essentially unchanged, 
indicates qualitatively that also for twisted Wilson fermions with small /i an 
Aoki phase has to be expected (compare [Hj). 

It is interesting to analyze the case of a twisted chiral operator, i.e., we 
now demand that the lattice Dirac operator Dq obeys the Ginsparg-Wilson 
equation 

75D0 Do 75 = aDo 75 iAo. ( 7 ) 

Since a chiral Dq has exact zero modes with definite chirality = 

it is more suitable for analyzing the interplay between the twisted 
mass and topology. A twisted chiral Dirac operator provides a clean setting 
for such an analysis and the twisted Wilson Dirac operator is expected to 
approximate the chiral behavior for sufficiently smooth gauge helds. 

It is straightforward, that adding a twisted mass term to a chiral Dq turns 
" 0 ^ into an eigenmode of D{^) with eigenvalue A = ±f/i. However, more 
information can be extracted from the Ginsparg-Wilson equation. Using the 
dehnition one has Dq = D{fi) — Inserting this expression into 0 

one hnds after a few lines of algebra (use (0 and D{—fi) = D{fi) — z2^75), 

= aD{fi)D{^)^ — afj? . 

Multiplying this equation with an eigenvector V’A from the right and with 
"01 from the left, it turns into an equation for the eigenvalue A, 

A + A* = aXX* - . (8) 
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Figure 2; Graphical illustration of our results for spectrum (l.h.s. plot) and 
75 matrix element (r.h.s.) for twisted chiral fermions. 


Setting \ = X + iy one finds that this is the equation for a circle in the 
complex plane. This “stretched Ginsparg-Wilson circle” has center 1 /a 
and radius Thus, in addition to obeying the bound (jHl), the 

eigenvalues of a twisted chiral Dirac operator are restricted to the stretched 
Ginsparg-Wilson circle. 

Our results for the spectrum of twisted chiral fermions are represented 
in the l.h.s. plot of Fig. |21 The eigenvalues are restricted to the two arcs 
drawn as full curves. The zero-modes of the untwisted chiral operator are 
mapped to the edges ±i/x of the arcs (the doublers of the zero modes end 
up at the other edges at 2 /a ± iy). The index theorem for zero twist turns 
into a “remnant index theorem” for the edge modes, 

Q = N_ - N+, 


where N± is the number of eigenvalues A = ±i//. The corresponding eigen¬ 
vectors have chiralities ±1. 

The largest imaginary parts that can be reached by the eigenvalues re¬ 
stricted to the arcs are ImA = i-^/l/a^ -|- fj?. According to Q, this implies 
a lower bound for the 75 matrix element, and we conclude 




< KV^AiTsV^a)! < 1 ■ 


The upper bound is reached only by the topological edge-modes (and their 
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Figure 3: Properties of the topological edge mode as a function of We 
show the eigenvalue, the 75 matrix element and the inverse participation 
ratio IPR (left to right). The dotted lines are the results for a twisted chiral 
Dirac operator, the symbols are for twisted Wilson fermions. 


doubler partners). The relation Q between the chirality s-nd ImA 

is illustrated in the r.h.s. plot of Fig. |2J 

Analyzing the case of a twisted chiral fermion, we have seen that topology 
plays an important role for the spectrum also for finite /r. In particular the 
topological edge modes are protected by chiral symmetry and are unchanged 
when varying fi. However, Wilson fermions are not chiral and one expects 
that the connection between topology, low-lying spectrum and chirality of 
eigenmodes is only approximate. An instance of this approximate realization 
is apparent in the fact that the topological mode in the r.h.s. plot of Fig. Q 
is indeed the eigenvalue closest to the real axis. It is, however, not located 
at exactly —ifi, as would be the case for a twisted chiral Dirac operator. 

In order to study the behavior of the edge modes we computed the 
spectrum and eigenvectors of twisted Wilson fermions in the background 
of a discretized smooth instanton with radius /? = 4a on a 16^ lattice (see 
|14) for details of the discretization). In Fig. 0 we show, as a function of 
/r, the eigenvalue of the edge mode (l.h.s. plot), its chirality (center), and 
its inverse participation ratio (r.h.s.). The inverse participation ratio is 
dehned as is a measure for the localization of the 

eigenmode with large values corresponding to localized modes while small 
values indicate spread-out modes. In all three plots the symbols are the 
numerical data, while the dotted lines are the analytical results for twisted 
chiral fermions. The plots show a clear difference between the data and the 
expectation for chiral fermions. This indicates that twisted Wilson fermions 
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Q = -1, real Polyakov loop Q = +1, complex Polyakov loop 


^ = 0 ^ = 0.02 n = 0 n = 0.02 



Figure 4: Spectra at // = 0, 0.02 for configurations in the deconfined phase. 


do not give rise to a clear separation of topological modes and bulk modes, 
which is only restored in the limit /r ^ 0. Repeating the same analysis 
with the chirally improved Dirac operator we found that there the data 
points fall on the chiral lines, demonstrating that for this case the topological 
edge mode is essentially protected by chiral symmetry. We remark, that 
since the chirally improved Dirac operator is only an approximate solution 
of the Ginsparg-Wilson equation, in this case the protection of the edge 
mode is also only approximate. 

A final question which we address in this letter is the role of the twisted mass 
term for configurations in the deconfined phase. There the spectrum acquires 
a gap near the origin and, according to the Banks-Casher formula |lb| . the 
chiral condensate vanishes. In Fig. |3] we show the effect of the twist term 
by comparing spectra with fi = 0 and /i = 0.02 for quenched configurations 
in the deconfined phase (Liischer Weisz action, (3 = 8.45,20^ x 6). The 
l.h.s. pair of spectra is for a configuration with real Polyakov loop where 
a larger gap opens up dl, while the right-hand side is for a configuration 
with complex Polyakov loop. The plots show that in both cases, at least 
for moderately small the twisted mass term affects only the topological 
mode, while the bulk of the spectrum is essentially unchanged. 
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